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Abstract 

New geometric objects on null thin layers are introduced and their importance for crossing 

null-like shells are discussed. The Barrabes-Israel equations are represented in a new geometric 

form and they split into decoupled system of equations for two different geometric objects: 

^. . tensor density G a b and vector field I. Continuity properties of these objects through a crossing 

■^j- ' sphere are proved. In the case of spherical symmetry Dray-t'Hooft-Redmount formula results 

00 from continuity property of the corresponding object. 

o" 

O ; 1 Introduction 

Self gravitating matter shell (see |151 118] ') became an important laboratory for testing global 
properties of gravitational field interacting with matter. Models of a thin matter layer allow 
us to construct useful mini-superspace examples. Toy models of quantum gravity, started by Dirac 
5-H i [S] , may give us a deeper insight into a possible future shape of the quantum theory of gravity (see 

Pf" El)- Especially interesting are null-like shells, carrying a self-gravitating light-like matter (see 

|ll)MllllT21ll3p . Classical equations of motion of such a shell have been derived by Barrabes and 
Israel in their seminal paper [3J- Junction conditions for general hypersurfaces in spacetime are 
also given in |22j . 

A complete Lagrangian and Hamiltonian description of the theory of self-gravitating light-like 
matter shell, which is no longer spherically symmetric, was given (in terms of gauge-independent 
geometric quantities) in J7j . For this purpose the notion of an extrinsic curvature for a null- like hy- 
pcrsurfacc was discussed and the corresponding Gauss-Codazzi equations were proved. These equa- 
tions imply Bianchi identities for spacctimes with null-like, singular curvature. Energy-momentum 
tensor density of a light-like matter shell is unambiguously defined in terms of an invariant matter 
Lagrangian density. Noethcr identity and Bclinfantc Roscnfcld theorem for such a tensor density 
was also proved. Finally, the Hamiltonian dynamics of the interacting system: "gravity + matter" 
was derived from the total Lagrangian, the latter being an invariant scalar density. 

Starting from the action functional for a single spherical shell due to Louko, Whiting and 
Friedman |21j . Hajicek and Kouletsis generalized it for any number of spherically symmetric null 
shells, including the cases, when the shells intersect |12| . 

In this paper we consider a general non-symmetric case of two crossing null shells. It occurs 
that the geometric objects on the null shells arc continuous through an intersecting sphere due to 
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the observation that "jump of the jump" vanishes (see Lemma l4.ip . This implies that the dynamics 
of the crossing shells is described by the equations for a single shell plus continuity property across 
intersecting sphere. 

We also discuss a special case of spherical symmetry. In particular, we give a simple argument 
(in the case of spherical symmetry) for triviality of the whole "ADM-momentum" tensor density 
G a b which implies that the corresponding energy-momentum tensor density r a b of a light-like 
matter shell is vanishing. 

Geometry of a single shell introduced in ^7j is completed by an extra object — a null vector 
field /, which is always well defined on a null shell and does not vanish in the case of spherical 
symmetry. Roughly speaking, in the case of a null shell the "ADM-momentum" tensor density 
G afc (which is well defined for any non-degenerate surface S) splits into two geometric objects: 
a tensor density G% and a null vector I a . They contain a similar information as the jump of a 
"transverse" extrinsic curvature K. a b in Barrabes-Israel approach. 

The dynamical system constituted of two spherically symmetric null shells has been studied 
in [IH- The shells at intersection sphere S x exchange energy according to the Dray-t'Hooft- 
Rcdmount formula jH] I25| . We show that the continuity of the metric (around intersection sphere) 
implies the continuity of the vector field I through S x on both shells. Moreover, in the case of 
spherical symmetry we show that the continuity of / gives the Dray-t'Hooft-Redmount formula. 
This means that our new object should be useful in generalizations of the Dray-t'Hooft-Redmount 
formula for the case of crossing two null shells without any symmetry. 

2 Geometry of a single null shell 

2.1 Geometry of a null hypersurface and Gauss— Codazzi constraints 

A null hypersurface in a Lorcntzian spacctime M is a three-dimensional submanifold S C M such 
that the restriction g a b of the spacctime metric g^ to S is degenerate. 

We shall often use adapted coordinates, where coordinate a; 3 is constant on S. Space coordinates 
will be labeled by k, I — 1, 2, 3; coordinates on S will be labeled by a, b = 0, 1, 2; finally, coordinates 
onSj := Vt n S (where V t is a Cauchy surface corresponding to constant value of the "time- like" 
coordinate a; = t) will be labeled by A, B — 1, 2. Spacctime coordinates will be labeled by Greek 
characters a, (3, ^i, v. 

The non-degeneracy of the spacctime metric implies that the metric g a b induced on S from the 
spacetime metric g^ v has signature (0,+,+). This means that there is a non- vanishing null-like 
vector field X a on S, such that its four-dimensional embedding X^ to M (in adapted coordinates 
X 3 = 0) is orthogonal to S. Hence, the covector X v = X^g^ u = X a g ai/ vanishes on vectors tangent 
to S and, therefore, the following identity holds: 

X a g ab = . (1) 

It is easy to prove (cf. ^HI) that integral curves of X a , after a suitable reparameterization, are 
geodesic curves of the spacetime metric g^ v . Moreover, any null hypersurface S may always be 
embedded in a one-parameter congruence of null hypersurfaces. 

We assume that topologically we have S = IR 1 xS 2 . Since our considerations are purely 
local, we fix the orientation of the M 1 component and assume that null-like vectors X describing 
degeneracy of the metric g a b of S will be always compatible with this orientation. Moreover, we shall 



always use coordinates such that the coordinate x° increases in the direction of X, i.e., inequality 
X(x°) = X° > holds. In these coordinates degeneracy fields are of the form X = /(d — n A dA), 

~AB 

where / > 0, riA = go A aud we rise indices with the help of the two-dimensional matrix g . 
inverse to gAB- 

If by A we denote the two-dimensional volume form on each surface x° = const.: 



A := yMet g AB , (2) 

then for any degeneracy field X of g a b the following object 

A 



I'X 



X(x°) 



is a well defined scalar density on S according to ^7]. This means that vx := vxdx A dx 1 A dx 2 
is a coordinate-independent differential three-form on S. However, vx depends upon the choice of 
the field X. 

It follows immediately from the above definition that the following object: 

A = v x X 

is a well defined (i.e., coordinate-independent) vector density on S. Obviously, it does not depend 
upon any choice of the field X : 

A = \(d - n A d A ) . (3) 

Hence, it is an intrinsic property of the internal geometry g a b of 5*. The same is true for the diver- 
gence d a A a , which is, therefore, an invariant, AT-independent, scalar density on S. Mathematically 
(in terms of differential forms) , the quantity A represents the two- form: 

L := A a (d a J dx° A dx 1 A dx 2 ) , 

whereas the divergence represents its exterior derivative (a three-from): dL := (<9 a A a ) dx° A dx 1 A 
dx 2 . In particular, a null surface with vanishing dL is called a non-expanding horizon (see |2]). 

Both objects L and vx may be defined geometrically, without any use of coordinates. For 
this purpose we note that at each point x £ S, the tangent space T X S may be quotiented with 
respect to the degeneracy subspace spanned by X. The quotient space carries a non-degenerate 
Riemannian metric and, therefore, is equipped with a volume form u> (its coordinate expression 
would be: u> — A dx 1 /\dx 2 ). The two- form L is equal to the pull-back of to from the quotient space 
to T X S. The three-form vx may be defined as a product: vx =«AL, where a is any one-form 
on S, such that < X, a >= 1. 

The degenerate metric g ao on S does not allow to define via the compatibility condition \7g = 0, 
any natural connection, which could be applied to generic tensor fields on S. Nevertheless, there 
is one exception: it was shown in |17| that the degenerate metric defines uniquely a certain covari- 
ant, first order differential operator. The operator may be applied only to mixed (contravariant- 
covariant) tensor density fields H a b, satisfying the following algebraic identities: 

n a b x b = o , (4) 

Hab = H-ba , (5) 



where H a b '■= 9acH. c b- Its definition cannot be extended to other tensorial fields on 5". Fortunately, 
the extrinsic curvature of a null-like surface and the energy-momentum tensor of a null-like shell 
are described by tensor densities of this type. 

The operator, which we denote by V a , is defined by means of the four-dimensional metric 
connection in the ambient spacetime M in the following way: Given H a b , take any its extension 
H MI/ to a four-dimensional, symmetric tensor density, "orthogonal" to S, i.e. satisfying H ±iy = 
("_L" denotes the component transversal to S). Define V a H a 6 as the restriction to S of the four- 
dimensional covariant divergence V^H^. It was shown in JJj that ambiguities which arise when 
extending three-dimensional object H° b living on 5* to the four-dimensional one, cancel finally and 
the result is unambiguously defined as a covector density on S. It turns out, however, that this 
result does not depend upon the spacetime geometry and may be defined intrinsically on S as 
follows: 

V a H a 6 = <9 a H a b — -H ac g ac .b , 

where g a c,b '■— dbg a c, a tensor density H a 6 satisfies identities (@J and (O, and moreover, H ac is any 
symmetric tensor density, which reproduces H a h when lowering an index: 

U a b = H ac g cb . (6) 

It is easily seen, that such a tensor density always exists due to identities 10} and JHJ), but the 
reconstruction of H ac from H a 6 is not unique, because H ac + CX a X c also satisfies (J5J) if H ac does. 
Conversely, two such symmetric tensors H ac satisfying © may differ only by CX a X c . Fortunately, 
this non- uniqueness does not influence the value of ©■ Hence, the following definition makes sense: 

V a H a 6 := d a n a b - ^U ac g ac , b • (7) 

The right-hand-side does not depend upon any choice of coordinates (i.e., transforms like a genuine 
covector density under change of coordinates). 

To express directly the result in terms of the original tensor density H a fc , we observe that it has 
five independent components and may be uniquely reconstructed from H° A (2 independent com- 
ponents) and the symmetric two-dimensional matrix Hab (3 independent components). Indeed, 
identities (@J and |[5j may be rewritten as follows: 

H A B = ~g AC n CB -n A K° B , (8) 

H° = VL\n A , (9) 

n B =(g_ BC n CA -n B n\)n A . (10) 

The correspondence between H a b and (H%,H^b) is one-to-one. 

To reconstruct H ab from H a b up to an arbitrary additive term CX a X b , take the following, 
coordinate dependent, symmetric quantity: 

F AB := g AC n CD9 DB - n A H° c f B n B H° c f A , (11) 



F 0A 



= H U C 3 =:F AU , (12) 

= 0. (13) 



It is easy to observe that any H ah satisfying © must be of the form: 

n ab =F ab + H 0Q x a x b _ ( U ) 

The non-uniqueness in the reconstruction of H is, therefore, completely described by the arbi- 
trariness in the choice of the value of H 00 . Using these results, we finally obtain: 

V a H a b := <9 a H a b — —H ac gac,b = d a H a b — -^F ac g ac ^ 

= d a II a b - \ (2H° A n\ - H AC ~g AC b ) . (15) 

The operator on the right-hand-side of i|15fl is called the (three-dimensional) covariant derivative 
of H a 6 on S with respect to its degenerate metric g a b- It was proved in J7| that it is well defined 
(i.e., coordinate- independent) for a tensor density H",, fulfilling conditions Q and (J5J. It was 
also shown that the above definition coincides with the one given in terms of the four-dimensional 
metric connection and due to 10, it equals: 

V^H* = d^\ - \n^g, lX . b = d a H a b - ^U ac g ac , b , (16) 

and, whence, coincides with V a H a 6 defined intrinsically on S. 

To describe exterior geometry of S we begin with covariant derivatives along S of the "orthog- 
onal vector X". Consider the tensor S7 a X^. Unlike in the non-degenerate case, there is no unique 
"normalization" of X and, therefore, such an object does depend upon a choice of the field X . 
The length of X vanishes. Hence, the tensor is again orthogonal to S, i.e., the components corre- 
sponding to [i — 3 vanish identically in adapted coordinates. This means that V a X b is a purely 
three-dimensional tensor living on S. For our purposes it is useful to use the "ADM-momentum" 
version of this object, defined in the following way: 

Q\(X) := -s {vx (V b X a - 5tV c X c ) + S^d c A c } , (17) 

where s := sgn<? 03 = ±1. Due to above convention, the object Q a b (X) feels only external orienta- 
tion of S and does not feel any internal orientation of the field X. 

Remark: If S is a non- expanding horizon, the last term in the above definition vanishes. 

The last term in i|17f) is -^-independent. It has been introduced in order to correct alge- 
braic properties of the quantity vx (V^X" — S b V C X C ): it was shown in [T7I that Q a b satisfies 
identities (@J-(0) and, therefore, its covariant divergence with respect to the degenerate metric 
g a b on S is uniquely defined. This divergence enters into the Gauss-Codazzi equations, which 
relate the divergence of Q with the transversal component Q\ of the Einstein tensor density 
Q^ = y/\ dct g\ {R% — S^R). The transversal component of such a tensor density is a well de- 
fined three-dimensional object living on S. In coordinate system adapted to 5 1 , i.e., such that the 
coordinate x 3 is constant on S, we have Q\ = Q\. Due to the fact that Q is a tensor density, 
components Q\ do not change with changes of the coordinate x 3 , provided it remains constant on 
S. These components describe, therefore, an intrinsic covector density living on S. 

Proposition 1. The following null-like-surface version of the Gauss-Codazzi equation is true: 
V a Q a b (X) + sv x d b {^f\ = -G\ ■ (18) 



We remind the reader that the ratio between two scalar densities: d c A c and vx, is a scalar 
function. Its gradient is a covector field. Finally, multiplied by the density Vx, it produces an 
intrinsic covector density on S. This proves that also the left-hand-side is a well defined geometric 
object living on S. The equation l|18|) is closely related to Raychaudhuri [21] equation for the 
congruence of null geodesies generated by the vector field X. 

2.2 Bianchi identities for spacetimes with distribution valued curvature 

In this paper we consider a spacetime M with distribution valued curvature tensor in the sense of 
Taub [27j . This means that the metric tensor, although continuous, is not necessarily C^-smooth 
across S: we assume that the connection coefficients T* may have only step discontinuities (jumps) 
across S. Formally, we may calculate the Ricmann curvature tensor of such a spacetime, but 
derivatives of these discontinuities with respect to the variable x 3 produce a (5-like, singular part 
of R: 

sing(#)V« - (WJ - S 3 K [T^)) d(x 3 ) , (19) 

where by S we denote the Dirac distribution (in order to distinguish it from the Kronecker symbol 
<5) and by [/] we denote the jump of a discontinuous quantity / between the two sides of S. The 
above formula is invariant under smooth transformations of coordinates. There is, however, no sense 
to impose such a smoothness across S. In fact, the smoothness of spacetime is an independent 
condition on both sides of S. The only reasonable assumption imposed on the differentiable 
structure of M is that the metric tensor — which is smooth separately on both sides of S — remains 
continuous across S. Admitting coordinate transformations preserving the above condition, we 
loose a part of information contained in quantity 1(1911 . which becomes now coordinate-dependent. 
It turns out, however, that another part, namely the Einstein tensor density calculated from 
(|19fl . preserves its geometric, intrinsic (i.e., coordinate- independent) meaning. In case of a non- 
degenerate geometry of S, the following formula was used by many authors (see |15l 1181 171 151 15] 1: 

sing(S)^ = G^S(x 3 ) , (20) 

where the "transversal-to-S"' part of C" vanishes identically: 

G ±,y = 0, (21) 

and the "tangent-to-S" part G ab equals to the jump of the ADM-momentum 1 Q ab of S between 
the two sides of the surface: 

G ab = { Q ab} ( 22 ) 

This quantity is a purely three-dimensional, symmetric tensor density living on S. When multiplied 
by the one- dimensional density S(x 3 ) in the transversal direction, it produces the four- dimensional 
tensor density Q according to formula 1)20(1. 

In the case of our degenerate surface S it was shown in ^JJ that formulae 1(20(1 and (|21|l remain 
valid also in this case. In particular, the latter formula means that the four-dimensional quantity 

1 Q a = yj | det g c d | (g ab trK — K ab ), where g ab is the inverse three-metric and K ab is an extrinsic curvature. 



QP" reduces in fact to an intrinsic, three-dimensional quantity living on S. However, formula (|22[) 
cannot be true, because — as we have seen — there is no way to define uniquely the object Q ab 
for the degenerate metric on S. Instead, we are able to prove the following formula: 

G\ = [Q\{X)\ , (23) 

where the bracket denotes the jump of Q a b (X) between the two sides of the singular surface. This 
quantity does not depend upon any choice of X and the singular part sing(C?) a & of the Einstein 
tensor is well defined. We will show in the sequel that the missing component G 00 can be recovered 
in another geometric object, which is presented in the next Section. 

Remark: Otherwise as in the non-degenerate case, the contravariant components G ab in formula 
(|20|) do not transform as a tensor density on S. Hence, the quantity defined by these components 
would be coordinate-dependent. According to J22J, G becomes an intrinsic three-dimensional 
tensor density on S only after lowering an index, i.e., in the version of G a &. This proves that C" 
may be reconstructed from G a b up to an additive term CX^X" only. We stress that the dynamics 
of the shell is unambiguously expressed in terms of the gauge-invariant, intrinsic quantity G a &. 

We conclude that the total Einstein tensor of our spacetime is a sum of the regular part 2 reg(£7) 
and the above singular part sing(CJ) living on the singularity surface S. Thus 

g^=reg(gr u + smg(gr i/ , (24) 

and the singular part is given up to an additive term CX fJl X L ,S(x 3 ). The following four- dimensional 
covariant divergence is unambiguously defined: 

= V„S" C = d^ c - g^ a T» c = d^ c - \g» x g^, c . (25) 

It is proved in |17| that this quantity vanishes identically and the total singular part of the Bianchi 
identities reads: 

sing (V^J = ([reg(0) X c ] + V a G a 6 ) S(x 3 ) = , (26) 

and vanishes identically due to the Gauss-Codazzi equation <|18|) . when we calculate its jump 
across S. Hence, the Bianchi identity V^g^ = holds universally (in the sense of distributions) 
for spacetimes with singular, light-like curvature. 

It is worthwhile to notice that the last term in definition (|17f) of the tensor density Q of S is 
identical on its both sides. Hence, its jump across S vanishes identically. This way the singular 
part of the Einstein tensor density $Zty reduces to: 

G\ = [Q\] = svx ({V b X a } - 6Z{V C X C }) . (27) 

2.3 Energy-momentum tensor of a light-like matter. Belinfante-Rosenfeld 
identity 

The interaction between a thin light-like matter-shell and the gravitational field is described in 
|17j . In particular, all the properties of such a matter are derived from its Lagrangian density L, 

2 The regular part is a smooth tensor density on both sides of the surface S (calculated for the metric g separately) 
with possible step discontinuity across S. 



which depends upon (non-specified) matter fields z K living on a null-like surface S, together with 
their first derivatives z K a := d a z K and — of course — the (degenerate) metric tensor g ab of S: 

L = L(z K ;z K a ;g ab ) . (28) 

We assume that L is an invariant scalar density on S. Similarly as in the standard case of canonical 
field theory, invariance of the Lagrangian with respect to reparameterizations of S implies impor- 
tant properties of the theory: the Bclinfantc-Rosenfeld identity and the Noether theorem, which 
will be discussed in this Section. To get rid of some technicalities, we assume in this paper that the 
matter fields z K are "spacetime scalars", like, e.g., material variables of any thermo-mechanical 
theory of continuous media (see, e.g., [HIGHI)- This means that the Lie derivative Lyz of these 
fields with respect to a vector field FonS coincides with the partial derivative: 

(C Y z) K = z K a Y a . 

The following Lemma characterizes Lagrangians which fulfill the invariance condition: 

Lemma 2.1. Lagrangian density \2tty concentrated on a null hypersurface S is invariant if and 
only if it is of the form: 

L = v x f(z;C x z;g), (29) 

where X is any degeneracy field of the metric g a b on S and /(• ; ■ ; ■) is a scalar function, homo- 
geneous of degree 1 with respect to its second variable. 

Remark: Because of the homogeneity of / with respect to Cxz, the above quantity does not 
depend upon a choice of the degeneracy field X. 

Dynamical properties of such a matter are described by its canonical energy-momentum tensor 
density, defined in a standard way: 

ell 
T\ := ^ ir z K b - 5%L . (30) 

®z a 

It is "symmetric" in the following sense: 

Proposition 2. Canonical energy-momentum tensor density T a b constructed from an invariant 
Lagrangian density fulfills identities O) and |3J), i.e., the following holds: 

T\X b = and T ab = T ba . (31) 

In case of a non-degenerate geometry of S, one considers also the "symmetric energy-momentum 
tensor density" r ab , defined as follows: 

T ab ._ 2 ^±L (32) 

og a b 

In our case the degenerate metric fulfills the constraint: det g ab = 0. Hence, the above quantity 
is not uniquely defined. However, we may define it, but only up to an additive term equal to the 
annihilator of this constraint. It is easy to see that the annihilator is of the form CX a X b . Hence, 
the ambiguity in the definition of the symmetric energy-momentum tensor is precisely equal to 



the ambiguity in the definition of T ab , if we want to reconstruct it from the well defined object 
T a b. This ambiguity is cancelled, when we lower an index. The next theorem says that for field 
configurations satisfying field equations, both the canonical and the symmetric tensors coincide 3 . 
This is an analog of the standard Belinfante-Roscnfcld identity (see 0]). Moreover, Noethcr 
theorem (vanishing of the divergence of T) is true. We summarize these facts in the following: 

Proposition 3. If L is an invariant Lagrangian and if the field configuration z K satisfies Euler- 
Lagrange equations derived from L: 

dL n dL 

d^~ d «d^-r°> (33) 

then the following statements are true: 

1. Belinj "ante- Rosen] "eld identity: canonical energy-momentum tensor T a \, coincides with (minus 

- because of the convention used) symmetric energy-momentum tensor t : 

T\ = -r ac g cb , (34) 

2. Noether Theorem: 

V a T\ = . (35) 

It is shown in |17j that the Einstein equations for the singular part: 

G\ = 8nT a b (36) 

can be derived from an action principle and they contain an intrinsic part of the Barrabes-Israel 
equations in mixed (contravariant-covariant) tensor density representation. Let us notice that if we 
assume vacuum Einstein equations outside surface S then, in particular, they imply reg(t/) J c = 
which gives compatibility of Q26(l with l|35|l . 

Remark: We may also include a regular matter part into the action and we obtain that the regular 
part of the energy momentum tensor density is no longer vanishing. In that case our null singular 
matter fulfills the following equation: 

sing (V M T£) = ([regCO-M + V a r a b ) S(x 3 ) = , (37) 

where T^ is the symmetric energy-momentum tensor density of the whole matter surrounding our 
shell S. If reg(T) A " / is derived form the (regular part of) Lagrangian then cq. I|36|l may be also 
considered as a generalized Noether theorem for the the full (regular + singular) Lagrangian of 
matter. 



3 In our convention, the energy is described by formula: H = T% = Pk°z k — L > 0, analogous to H = pq — L 
in mechanics and well adapted for Hamiltonian purposes. This convention differs from the one used in |23l . where 
the energy is given by Too- To keep standard conventions for Einstein equations, we take standard definition of the 
symmetric energy-momentum tensor r a (,. This is why Bclinfantc-Rosenfcld theorem takes form r a t, = —T a ^. 



3 Canonical null vector on a single shell 

Let us rewrite the Ricci tensor: 

Rnu = d\Y^ u - <9( M r„ )A + r CTA r£„ - r M(T r^ A , (38) 

in terms of the following combinations of Christoffel symbols: 

A%, := T^ - 5^r: )K . (39) 

We have: 

R lw = d x A^-A^A: x + ^ x A: a . (40) 

The terms quadratic in A's may have only step-like discontinuities. The derivatives along 5* are 
thus bounded and belong to the regular part of the Ricci tensor. The singular part of the Ricci 
tensor is obtained from the transversal derivatives only. In our adapted coordinate system, where 
x 3 is constant on S, we obtain: 

sing(iV) = daA%, = 6(a?)[A%,} , (41) 

where by 8 we denote the Dirac delta-distribution and by square brackets we denote the jump of 
the value of the corresponding expression between the two sides of S. Consequently, the singular 
part of Einstein tensor density reads: 

sing(0"„) := \/M sing (V„ - ±r\ = 6(x 3 )G» v , (42) 

where 

G^ := y/\g\ {staT ^ W 9 ) [AU = [QM , (43) 

Q^ ■= V\g\ {gTsf* - \iTtf*) A% , (44) 

and explicit formulae for Q^ „ are given in Appendix [BJ It was also shown in 17 that the con- 
travariant version of this quantity: 

smg^r = iQns(x 3 ) , 

is coordinate-dependent and, therefore, does not define any geometric object. Let us observe that 
G ab := [Q ab ] is not well defined intrinsic tensor density on S in contrast to G a b = [Q a &], as was 
shown in Appendix A of |17j . However, one can extract the following object: 

pOO 

7 ° :=sX °x^' (45) 

which is well defined because of the following 

10 



Proposition 4. The vector field I defined by \45\j does not depend on the choice of the field X 
and coordinate x° , hence it is a well defined intrinsic object on the null surface S. 

Proof. Let us express the component G 00 in terms of the objects which arise in (l+2+l)-decomposition 
of spacetime (see Appendix [5J): 

G 00 =[Q 00 H.9 03 [0°3]+.9 06 [Q 



= ^{-\d,\n\]+m\w b ])-s{^X b + ±n 

= - jj[d 3 \] = -sY»{d,\] , (46) 

where the last equality holds because tangent to S derivatives d a X are continuous, hence [<9 a A] = 0. 
The transformation laws, introduced in |16| and given in Appendix^ imply that 

pOO 

sXa XOAo=- Y ^ lnX ^ Xa 

is not dependent on the choice of the basis X, 8a, Y at the point ieS. More precisely, for any 
two tetrads X, d Al Y and X,d^,Y related by (E2JHIE3, J75J) we get [Y(ln\)]X = [Y(lnX)]X. □ 

We also have G^Y^Y" = G 00 because G^X^ = (cf. J2TJ and Appendix 01 . 

Remark: One can define a symmetric tensor density W := I ® A — A <g> / on S. However, there is 
no possibility to include object W ab into G ab unless G a ;, = 0. Moreover, if G a b is vanishing (which 
happens for spherical symmetry cf. Prop. EJ), one can check from Bianchi identities V ^Q^v = 
that 

for any extension I 11 which is tangent to S. Unfortunately, this equation is not intrinsic on 5*. 

The equation l|36|l cannot be completed by the equality G 00 = 8ttt 00 on the tensor density 
level because nor G ab neither r ab are geometric objects 4 on S. On the other hand, the definition 
(I45|l allows to complete singular Einstein equations 13(j|) in the following form: 

where the vector field P a defined as follows: 

pa '- =sXa xJk» (48) 

contains missing information about singular energy- momentum tensor density t^ v . 

Let us finish this section with the following observation: for non-degenerate surface S the tensor 
density G ab (given by (|22|l 'l is well defined. For the null shell S it splits into two objects: the tensor 
density G a £, defined by l|23|l and the null vector I a given by l|45|) . This means that the information 
about the jump of a "transverse" extrinsic curvature K, a b (in Barrabes-Israel approach) is contained 
in two different geometric objects - G a b and I a . 

4 In non-degenerate case both tensor densities are well defined. 

If 



4 Crossing shells 

Let us consider two shells intersecting each other along surface S x which is a sphere. One can 
imagine this situation with the help of Fig. ^ where one spherical coordinate is suppressed and 




Figure 1: Crossing shells 

the spheres are drawn as one-dimensional circles. 

Let us introduce a local coordinate system (v t x ,u) around S x , such that N u := {v, = uo} is 
the first shell and N v := {v = vq} is a second one. Hence S x = N u nN v . The metric takes the form 
similar to l|77J) but now both (transversal to S x ) coordinates u and v are null, i.e. corresponding 
level three-surfaces are degenerate. More precisely, 



9tiv 



n riA 


n A 


sM + m it>a 


n A 


9AB 


m A 


sM + m A nA 


m A 


ra A m,A 



(49) 
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which gives -J\ det g^\ = AM, and the contravariant four-metric takes the form 






m 

* M 


s 
A/ 


s M 


~~ AB . „n A m B +m A n B 


* M 


s 
M 


a M 






(50) 



where M > 0, s := sgn 5™ = ±1, gAB is the induced two-metric on surfaces {u = const, v = const} 

and g is its inverse (contravariant) metric. Both g and gAB are used to rise and lower indices 
A, B = 1,2 of the two- vectors n A and m A . 
Let us choose the null vector fields 



K 



d v — n 8a and L := d u — m ua 



which are tangent to N u or N v respectively and g{K, L) = sM. We can use the coordinates (v, x A ) 
on the first shell N u . On the second shell N v we have the coordinate system (x A , u). The canonical 
vector field / is well defined on both shells: 



I{K) - -§[£(lnA)]„ , I(L) = ~[KQn\)] v 



(51) 



where the index u or v corresponds to jump across first or second shell respectively. 

Several continuity properties of discontinuities across S x arc implied by the observation that 
jump of the jump vanishes which we explain below on the example of a real function of two variables. 
Let / be a function on an open set U C R 2 containing point (0, 0) such that / is smooth outside 
axes (corresponding to our crossing shells), i.e. / <G C k (U') for sufficiently large k > 2 and 

U'~U\ ({(x, y) G R 2 I x = 0} U {(x, y) G M 2 \ y = 0}) . 

Moreover, we assume that / is continuous across the axes with finite jumps of first normal deriva- 
tives. More precisely, the jump 



df 



d.r 



lim — (x,y) - lim -^-(x,y) 

x— >0+ OX x^Q- OX 



is well defined for 1/ ^ and splits into upper (positive y) and lower (negative y) parts. Under 
above assumptions we get 



Lemma 4.1. The jump 
'df 



lim 



d.r 



(y) 



lim 



is continuous across (0,0), i.e. 
df] 



O.r 



(y) 



and the similar property holds on x-axis. 
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Proof. Let us enumerate the quadrants of the plane: I, II, III, IV, i.e. I 



0,y > 0}, II — {(x,y) S 



< 0,2/ > 0}, III — > {(x,y) e 



-> {(a;, 2/) £ K 2 | x > 
< 0,2/ < 0}, IV — ► 



{(x, y) el 2 | x > 0,2/ < 0}, and the corresponding restrictions of the function / we denote by 
index e.g. the function / in the second quadrant we denote by f 11 . Continuity of / and its tangent 

derivatives across positive y-half-axis implies / 7 (0,2/) = f TI (0,y) and 4^ (0,2/) = 4) (0>2/) 
n = 1, 2, . . . k, where the boundary values of / and its derivatives are defined in an obvious way 
e.g. f I (Q,y) = lim 1E _ > o+ / 7 (x,2/). In particular, we have 

df 1 df H 

^- (o,y) = ^ (0,y) forj/>o, 

QfIV Q.III 

-± (0,y) = f- (0,y) iory<0. 
Passing to the limit at (0, 0), we get 



• ii 



• ii 



df 1 df 1 df 11 df 1 

%- (0,0):= lim^ (0,2/)= lim ^- (0,y =: ^ (0,0) 

dy y^o+ dy y->o+ dy dy 



and similarly 



• III 



QfIV QrlV Q.III Qrl 

%- 0,0 := lim ?L (0,2/)= hm^ P,y)=:^ (0,0). 
dy y^o- dy y^o- dy dy 

Finally, from the last two equations we get 

Qrl QrlV Qrll Q.III 

i M ~i (0 ' 0) =di (0 ' 0) -di (0 ' 0) 



which implies continuity of jump 



df 



Ox 



We can denote symbolically the result as 



across y = 
df 



a 



d.r 



0, i.e. jump of the jump at the crossing 



point vanishes. 

Using Lemma |4. II one can show the following 



Theorem 1. The continuity of the metric across null shells implies that the vector fields I{K) 
and I(L) are continuous across S x ■ 

Moreover, from Lemma PTT1 we get that G a b{K) on N u and G a b{L) on N v are also continuous 5 
across S x . 

Proof. From definition l|51|l of the null field K and l|49|) we have 



5 Although G a b (K) does not depend on the choice of the null field K, we keep this argument to distinguish the 
shells. Moreover, we should remember that the coordinates x a depend on the shell, i.e. (x a ) = (v,x A ) for N u but 
(x a ) = (u,x A ) for N v . 
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hence we apply Lemma T4. II for the function A. More precisely, we take 

f{x, y) := X(u = x + u Q , v = y + v ,x A ) 

with fixed coordinates x , hence the point x — 0, y = corresponds to the fixed point on S x with 
coordinates x A . 

For G a biK ) we observe that 

G\iK) = S A a [w b ] = ^jK a K c [d u g cb ] u 

which is implied by l|87|) and (|90() - (|91|l . Moreover, from l|91|) we get [w v ] = n^uu], hence it is 
enough to consider 

[wa] = -^gAB[d u n B } u 

implied by (|90(l . and using Lemma f4. II for the function / := n B we obtain the result. □ 

The above Theorem and the considerations from Section 2 imply that the dynamics of crossing 
shells is described by equations $36)) and \37jj which hold on both shells plus continuity property 
across S x ■ 

4.1 Spherically symmetric shells 

Proposition 5. For spherically symmetric null shell the tensor density G a b is vanishing. 

This implies that the dynamics of the spherical shell is very simple, i.e. T a b = 0, hence eqs. 
(I36|l - H37ll arc trivially satisfied but vector field / is not vanishing as we show in the sequel. 

Proof. From (J23, {E3 and iJ^TJ we get 

G a 6 = [Q\] = A a [w b ] (52) 

but spherical symmetry gives [wa] = and, moreover, l|91|l implies [wq] =0. D 

Let us check the value of / for the spherical null shell which arises from matching two Schwarzschild 
metrics along spherically symmetric null surface. 

gi = -ll — j Au 2 - 2dudn + rfdn , i = l,2, (53) 

where 

dfl := d6 2 + sin 2 8dtp 2 . 

We take u > for gi and u < for g%, and r^R, u) := R + ^m. This implies that the full metric 
is continuous in coordinates (u,i?) across the shell u = 0. More precisely, 

2m i \ , 2 r. , / , ^ , m i 



ffi|„=n = " ( 1 j du ~2du[dR+— -du) +R z dfl= -du 2 - 2dudR + R 2 dfl 
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1 _ ^l\ du 2 _ 2du (dR + ^-du) + R 2 dn = g 2 \ _ n . 
R J V R / u 



Moreover, if we choose null field X = JL then the transversal field may be chosen as Y — — -8- 
and A = rf sin. 9, hence 

y (lnA )| iOUn+U^i ** 

v 'l«=o du \ R J lu =° R 2 

and finally 

l-^g^X. (54) 

Next, for crossing two spherical null shells we may check the Dray-t'Hooft-Redmount formula [Bj. 
|25j as follows: firstly we apply Theorem^ which from continuity of the metric implies continuity 
of the vector field /, secondly we check that the vector field / is continuous through the crossing 
sphere iff the Dray-t'Hooft-Redmount formula is true. 

Theorem 2. // the shells are spherically symmetric than continuity of the vector field I gives 
Dray-t'Hooft-Redmount formula J 671) . 

Proof. Let us consider the full description of crossing spherically symmetric null shells which can 
be nicely given in Kruskal-Szekeres coordinates (instead of Eddington-Finkclstcin used in (JSSJ)- 
We assume four domains (cf. Fig. |2J equipped with the Schwarzschild metrics 

32m 3 ( r \ 

gi = — ?-exp -- — jduidVi + rt dO , i = 1,2,3,4, (55) 

r% \ Znii J 

where r{ = 2min(—UiVi) and the Kruskal function n is defined by its inverse k~ 1 (x) = (x — l)e x 
on the interval (0, oo) C R. One can easily check the following identity for the first derivative of 
k: 

k> = *! . (56) 

The four domains Mj (i = 1, 2, 3, 4) are matched together along null surfaces {x £ Mi \ Ui = a.;} C 
Mi and {x £ Mi \ Vi = Pi} C Mi, as is shown on Fig. [3 

The coordinates V\, V4 on the shell A^ 14 do not match but 

r = 2?7Jik(— a\Vi) = Im^K^—u^v^ (57) 

is the same on both sides and can be chosen as a coordinate on the surface -/V14. This equality 
also means that A is continuous across this shell. On the other hand the continuity of the term 
—^- exp ( —j^- ) duidvi across 7V 14 implies 



,:s 



exp I ) dwidf] 

r V 2m i 



— exp I I du4d«4 

r V 2to 4 
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f— 

second shell iV 12 



first shell N 23 




first shell N-44 



u 2 > a 2 , v 2 < 02 2 



U4 = ay 



second shell ^34 



Figure 2: Matching domains 



hence using (|56|l and l|57|) we obtain the transformation law between first derivatives of coordinates 
U\ and 114: 



d,U4 1 mi 

— — = — I exp 

du\ \ 777,4/ V 2toi 



I r \ dv\ m\a.4 
exp 



2m4j dvi 7714121 
Moreover, the null vector field X tangent to the first shell can be represented in M\ as follows 

d ( dr \~ d 

dr \ dv\ ) dvi 

and using 1)56(1 we have 

dv 1 

— — = — 2miai«'(— aiVi) = — 2m\ai («;(— aiui)exp[«;(— ai^i)]) , 
dv\ 



(58) 



hence 



X 



k{— ai^i) exp(/t(— a\V\)) d 



2a.\m\ dvi 

The transversal vector field 

Y = ai d 

fulfills normalization condition gi(X, Y) = 1. Moreover, using equality 

ai 9ri ri-2mi 

Y (ha A) = — — = 5 

2ra\r\ ou\ r\ 

and the similar one in A/4 we can check the formula ()54|l in new coordinate representation 



7 14 = -[y(lnA)]X = 



2rni — 2to4 



X. 
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Similar considerations for the first shell N23 give the following expression for the vector field H45|) : 
2m? — 2m 3 



'23 



-X. 



where now r = 2m2/t(— 012^2) = 2?ti3k(— 03^3) and 

d K(-a 2 V2)exp(K(-a 2 V2)) d K(-a 3 v 3 )exp(n(-a 3 v 3 )) d 



X 



dr 2a2Tn,2 dv 2 2a 3 m 3 dv 3 

We can compare 7x4 with ^23 across S'x by using the transformation law (cf. I|58|0 between V4 and 

"3 

dvj _ f3 4 m 3 
dv 3 (3 3 m A 

which is implied by continuity of the metrics g 3 and (74 across second shell JV34 (v 3 = (3 3 and 
Vi = Pa). Finally, we obtain 

2(m 2 - m 3 ) K(-a 3 ^ 3 ) exp(K(-Q 3 /3 3 )) d 

J 23V^3 P37 9 n rx 

r- 2 20137713 OU3 

2(m 2 - m 3 ) 2^7 ex P( 2^7) ftw 3 d 
r 2 2a 3 m 3 fi 3 ra\ dv4 

(m 2 -TO 3 )^4exp(2^-) a 
27-03/337/147713 <9t>4 

and 

2(mi-m 4 ) 2^I ex P(2^l) 9 

/ 14 (i; 4 = p A ) = — , 

hence ^23 = lu on S x implies 

(mi -m 4 )exp(2^) _ (m 2 - m 3 )f3 4 exp(^) 
2ra4rri4 2ra 3 f3 3 m 3 

or 

T T 

(mi - m4)a 3 /3 3 m 3 exp(-- ) = (m 2 - m 3 )a4/34m 4 exp(-- ) . (60) 

2777-3 2?Ti4 

Moreover, on S'x 

T T 

aiPi exp(-- — ) = 1 - /e(-a,-/?i) = 1 - - — 
Ami Zrrii 

which applied to (|60|) gives 

(mi — 7« 4 )(r — 2777.3) = (j n 2 ~ m 3)( r ~ 27714) , 
which is equivalent to Dray-t'Hooft-Rcdmount formula 

(7' — 2mi)(r — 27773) = (? ~ 2rri2)(r — 2777,4) ■ (61) 

D 

In the above proof we have restricted ourselves to the case of positive masses 777, and to the 
matching null surfaces which are not horizons. The analysis of possible special cases one can find 
in |11| but obviously the formula 1)61(1 remains valid for any special case. 
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A Transformation rules 

The triad {X, d A ) on S depends upon a particular (2 + l)-decomposition of S, given by the choice 
of the time coordinate x° on S. However, several objects constructed by means of the triad do 
not depend upon this choice and describe the geometry of S. To prove this independence, observe 
that we have the following transformation law: 



X = cX , (62) 

^B = C B 



B § = C£d A + f s X, (63) 



where (X,dg) is the new triad, corresponding to the new coordinate system (x a ) on S. The 
coefficient c may be obtained from the following equation: 



l = (d X °,X) (^dx A + ^dx°,cX) 

dx° A dSP 
dx A dx° 

hence, 



(dx° dx° 
\ dx° dx A 



n A 



(64) 



(65) 



On the other hand, we have: 

dx A dx° a\ „ 9a; 



-^n A 



hence, 



dx B dx B ) dx B 



dx A dx° 



d A + —^X , (66) 



B dx B dx B 

dx° 

f» - w- <68) 

The transformation law for qab'- 

9ab = Ca C b B 9 (9a + f A X, d B + f B X) = C^C^g AB (69) 

implies: 

A = AdetC/. (70) 
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In order to complete the triad (X, d A ) on 5 to a tetrad in M it is useful to choose a transverse 
field Y fulfilling the following "normalization conditions" : 



g(Y,x) 
g(X,d A ) 



o . 



(71) 
(72) 



These equations do not determine Y uniquely, but modulo an additive term proportional to X: a 
"gauge transformation" 



Y -^Y + hX , 



(73) 



with an arbitrary scalar field h is always possible. Extending coordinate x° from S to a neighbour- 
hood of S, we may choose the following transverse field: 



Y 



M & m " 



d A ) 



(74) 



We stress, however, that this particular choice of Y depends not only upon a (2 + l)-decomposition 
of S, but also on a (3 + l)-decomposition of M in a neighbourhood of S. Because of 11720 . the 
vectors X and Y span the bundle of vectors normal to S. 

The transformation law for Y, when passing from one to another (2 + l)-dccomposition of S, 
reads: 



Y = 



[Y - k A 3 A ) + hX , 



(75) 



where the scalar field h is arbitrary (it is determined by the extension of the (2 + l)-decomposition 
of S to a (3+ l)-decomposition of M), and the coefficients k are uniquely determined by equation 



fa = Cjf9Ack c 



(76) 



with /s given by (|68|l . Despite of the freedom in choice of Y, some geometric objects constructed 
with help of the tetrad (X, d A , Y) do not depend upon this choice and characterize only the 
geometry of S C M. 



B Structure of the singular Einstein tensor 

We are going to relate the coordinate-dependent quantity Q^ v with the external curvature Q a f, of 
S. We use the form of the metric introduced in |lfc>j : 



9^v 



n n A 


n A 


sM + m n A 


n A 


9AB 


m A 


sM + m A n A 


m A 


(f) 2 + m^m A . 



(77) 
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and 



-(^) 2 




n A „m A 
~W b M 


8 

M 


n A m A 

N 2 b M 


~AB 

9 


n A n B , „n A m B +m A n B 
N 2 ' '" M 


* M 


s 
M 




M 






(78) 



where M > 0, s := sgn<? 03 = ±1, gAB is the induced two-metric on surfaces {x° = const, 



zAB 



zAB 



const} and g is its inverse (contravariant) metric. Both g and gAB are used to rise and lower 



indices A, B = 1,2 of the two- vectors n A and m A . 



Formula (|77JI implies: y/\dctg^\ = AM. Moreover, the object A a defined by formula 0, takes 
the form A a = XX a . where A is given by formula J2J an d A" := do — n A dA- This means that we 
have chosen the following degeneracy field: A' 1 = (1, — n , 0). 

~AB 

For calculational purposes it is useful to rewrite the two-dimensional inverse metric g in 
three-dimensional notation, putting g := 0. This object satisfies the obvious identity: 

g ac g cb = S\ - X a S\ . 

Hence, the contravariant metric 1|78|) may be rewritten as follows: 



1 
N~2- 



ib _ ~ ab L -X a X b S ' 



M 



(m a X b + m b X a ) , 



(79) 



~aB 



where m a := g m,B, so that m := 0, and 



^3fi 



— X^ 

M 



It may be easily checked (see, e.g., [T^], P a g° 406) that covariant derivatives of the field X along 
S are equal to: 



zbc, 



V a X = -w a X - l ab g d c , 



where 



w a := -IT° , 



and 



Lb ■= -g(d b , V a X) = g{V a d b , X) = X^ ab 
Moreover, 



zab , 



d c A c = - Xg ab l ab = -Xfl ab = -XI , 

'1=9 lab- 

The following lemma was proved in |17j : 



(80) 

(81) 
(82) 

(83) 



where I = g l a b 
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Lemma B.l. The object Q a b is related to Q a b as follows: 

sQ\ = sQ\ - ^XIS\ + A a X b - S a b A c X c , (84) 

where \c ■= ^d c In f — 

Moreover, from definition (|i7|l and property (|8U[) one can check that 

sQ a b = X5 a b V c X c - XV b X a - 5 a b d c A c 

= -X5 a b (w c X c + l) + X(w b X a + ~g aC l cb ) + S a b M 

= \fX b + A a w b - 6 a b A c w c . (85) 

Remark: Formula l|85(l . together with l a bX b = = g ab X b , gives us the orthogonality condition 
Q a bX b = and symmetry of the tensor Q a b := g ac Q c b- 

Now, we would like to examine the properties of G^" = [Q^]. From continuity of the metric 
across S we obtain 

[l ab ] = sM[A 3 ab ] = sM[Tl b ] = X c [T cab ] = , (86) 

A a [A 3 b ] - S\A C [A 3 C ] 

A a {w b }-6 a b A c {w c }= s [Q a b }, (87) 



s[Q a b ] = A a [Al b ] - S\A C [A 3 C 



and 

[Q\] = (88) 

because sQ 3 3 = — \Xl and sQ 3 a = 0. 

Finally, the missing component [Q a 3] has the following form: 

[Q a 3 ] = sA a {-[8 3 In A] + m b [w b }} + MXf b [w b } . (89) 

We also have from 

K] = -X V 3 Fa ba ] = ^X b [g abt3 ] (90) 



that 



X a [w a ] = —[X a X b g ab , 3 ]=0. (91) 



Using these results from <|88|) one can easily check the property i|21[l 

G 33 = [Q 33 ] = g 33 {Q 3 3]+9 3b iQ\} = 0, 

G 3a = [Q 3a ] = g 33 [Q a 3 }+g 3b [Q a b ] = -^[X h Q a b } = 0, 

where we used the equality [Q a b] = [Q a b ] which is crucial to admit that the object G a b is a well 
defined geometric object on S. 
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C Gauss— Codazzi equations 

It was shown in |17j that 



sG\ = - sd b Q\ + ^sQ bc g bc , a + XdJ , (92) 



where we have used the formula 

sQ ab = \g g l cd + (A a g + A b g -g A c )w c . 
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